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DENSITY  CHANNEL  TRACKING  STUDIES  ON  PULSERAD 


I.  INTRODUCTION 

The  existence  of  a  mechanism  which  would  allow  electron  beams  to  track  preformed 
density  channels  is  a  key  element  for  long  range  propagation  of  Charged  Particle  Beams 
(CPBs).  Most  endoatmospheric  CPB  propagation  schemes  envision  multiple  beam  pulses 
boring  a  hole  in  the  atmosphere  as  the  beam  packet  propagates.  The  first  or  lead  pulse  of 
the  packet  propagates  in  a  uniform  density  gas  background.  As  it  propagates  it  deposits 
its  energy  in  the  gas,  heating  it  and  producing  a  reduced-density  channel.  If  succeeding 
pulses  follow  the  density  channel  produced  by  the  earlier  pulses  they  would  lose  less 
energy  in  the  lower  density  gas  than  they  would  outside  it.  Eventually  they  would  emerge 
from  the  end  of  the  preformed  channel,  become  lead  pulses  themselves  and  extend  the 
channel.  In  this  way  the  range  of  the  CPB  propagation  can  be  extended  beyond  the  uniform 
density  limit. 

An  experiment  using  a  single  pulse  accelerator  and  preformed  reduced-density  channels 
was  carried  out  to  test  the  basic  physics.  The  goal  of  the  research  was  to  experimentally 
test  a  tracking  mechanism  discovered  by  Welch  in  simulations  of  electron  beams  propa¬ 
gating  through  initially  un-ionized  air.[1]  A  summary  of  these  numerical  studies  and  the 
results  of  a  recent  tracking  experiment  performed  by  Bieniosek  at  McDonnell  Douglas 
Research  Labs  using  shallow  density  channels  is  given  in  reference  #1.  The  theory 
suggests  that  a  positive  tracking  force  exists  between  a  density  channel  and  an  electron 
beam  arising  from  asymmetric  plasma  current  flow  in  the  channel  induced  by  the  beam. 
The  strength  of  the  tracking  force  is  predicted  to  increase  with  the  beam  current  and  is  a 
function  of  the  channel  depth  and  the  beam  offset  relative  to  the  channel  axis.  The  IMeV 
beams  in  this  experiment  were  produced  by  the  1-MV  Pulserad  310  electron  beam 
generator[2],  and  the  reduced  density  channels  were  produced  by  the  technique  of  Laser 
Guided  Electric  Discharges  (LGED).[3]  The  electron  beam  pulses  were  about  35ns  long 
and  had  nominal  peak  beam  currents  of  9-1 2kA  in  the  first  experimental  run  (summer  1 987) 
and  currents  of  5-1  OkA  in  the  second  experimental  run  (summer  1988)  depending  on  the 
anode-cathode  (AK)  gap  selected  in  the  Pulserad  diode.  The  channels  were  offset  relative 
to  the  beam  axis  by  ±1  cm  in  the  first  run  and  by  Ocm,  ±1 .4cm  or  ±2.6cm  in  the  second  run. 
The  conducting  chamber  wall  radius  was  45cm  in  both  experiments  in  order  to  minimize 
wall  centering  effects. 

II.  DENSITY  CHANNEL  TRACKING  THEORY 

The  basis  of  the  density  tracking  effect  can  be  understood  from  an  analytic  theory 
developed  by  R.F.  Fen3ier.[4j  In  this  model  a  cylindrical  reduced-density  channel  is 
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surrounded  by  a  uniform,  higher  density  background  gas  distribution.  The  beam  is  initially 
displaced  off  the  channel  axis  by  a  distance  yc.  The  theory  is  based  on  three  key  elements: 

•  Ionization  of  the  background  gas  is  dominated  by  direct  beam  ionization.  This  yields 
an  ionization  fraction,  ne/ng  (ne=electron  density,  ng=local  gas  density),  on  the  beam 
timescale  that  is  independent  of  the  background  density  and  is  symmetric  about  the 
beam.  The  weak  plasma  fully  neutralizes  the  beam’s  space  charge. 

•  The  plasma  electrons  produced  by  the  beam  are  accelerated  and  heated  by  the  nearly 
uniform  axial  inductive  electric  field,  E,  produced  by  dleff/dt  where  leff  is  the  effective 
current  (beam  minus  plasma  current)  within  the  radius  of  the  beam.  The  plasma  electron 
temperature  in  this  weakly  ionized  gas  depends  on  the  reduced  field  strength  E/ng. 
Therefore  the  electron  temperature  in  the  lower  density  channel  will  be  higher  than  that 
in  the  surrounding  region. 

•  The  return  current  driven  by  the  inductive  electric  field  depends  on  the  electron-neutral 
collision  frequency  in  the  background  plasma.  For  electron  temperatures  up  to  several 
electron  volts  the  collision  frequency  increases  with  the  plasma  electron  temperature 
for  most  atmospheric  gases.  Thus  the  resistivity,  which  is  proportional  to  the  collision 
frequency,  will  be  higher  in  the  lower  density  central  region  than  in  the  higher  density 
regions  near  the  edge  of  the  channel.  The  return  current  driven  by  the  inductive  electric 
field  will  therefore  flow  preferentially  in  the  lower  resistivity  regions  near  the  edges  of 
the  channel.  When  the  beam  and  the  channel  are  offset  the  plasma  return  current 
distribution  will  be  asymmetric  about  the  channel  axis.  The  dipole  component  of  this 
plasma  current  interacting  with  the  beam  current  provides  the  centering  force  on  the 
beam  which  causes  the  beam  to  track  the  channel. 

Figure  1  shows  the  predictions  of  the  non-linear,  analytic  force  equation,  as  discussed  in 
Appendix  I.  For  a  beam  and  channel  similar  to  that  used  in  the  experiment  described 
below,  the  tracking  force  as  a  function  of  the  offset  between  the  channel  and  beam  is 
predicted  to  increase  from  zero  to  a  peak  near  the  channel  radius.  At  larger  offsets  the 
force  decreases.  The  tracking  length,  which  is  the  propagation  distance  necessary  for  the 
beam  to  cross  the  channel  axis,  is  independent  of  the  offset,  unless  the  offset  is  comparable 
to  or  larger  than  the  beam  or  channel  radii.  The  tracking  force  increases  linearly  with  the 
plasma  return  current  and  weakly  with  the  channel  gas  density.  Figure  2  shows  a  plot  of 
the  tracking  length  versus  relative  channel  gas  density.  The  tracking  length  monotonically 
increases  as  the  channel  depth  decreases.  The  beam  must  be  stable  over  the  tracking 
length  in  order  to  see  the  effect  clearly.  Tracking  lengths  for  the  pulserad  experiments 
range  from  under  40cm  to  over  1  m,  depending  on  the  beam  and  channel  parameters. 
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III.  EXPERIMENTAL  CONSIDERATIONS 


Several  effects  not  treated  in  the  theory  must  be  taken  into  account  for  a  viable  experiment 
is  to  be  performed.  The  beam  profile  is  degraded  by  scattering  off  the  background  gas, 
increasing  the  beam  radius.  The  beam  is  also  subject  to  the  hose  instability  which  can 
drive  the  beam  out  of  its  intended  trajectory.  Scale  lengths  for  these  effects  can  be 
calculated  and  compared  with  the  tracking  length. 

A  measure  of  the  effect  of  scattering  off  the  background  gas  is  the  Nordsieck  length,  Ln, 
which  measures  the  propagation  length  over  which  the  beam  radius  e-folds  due  to 
scattering.  A  rough  expression  for  the  Nordsieck  length  is: 


Ln  =  6.1 


(y  left)033 
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where  left  is  the  net  current  within  the  radius  of  the  beam,  i.e., 
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sion  for  the  Nordsieck  length  is  a  fit  to  numerical  simulations  by  Hughes  and  Godfrey.[5] 
Figure  3a  shows  the  Nordsieck  length  as  a  function  of  relative  channel  density  and  beam 
current.  The  beam  is  assumed  to  remain  inside  the  channel  which  significantly  increases 
Ln  for  a  low  density  channel.  The  tracking  length  for  even  a  low  current  beam  and  50% 
deep  channel  is  less  than  or  equal  to  the  Nordsieck  length.  Thus,  Nordsieck  expansion 
should  not  dominate  the  tracking  length  measurement. 

The  hose  instability  is  more  difficult  to  overcome  experimentally.  It  can  be  strong  enough 
to  drive  the  beam  off  its  initial  trajectory  and  overcome  a  positive  tracking  force.  Once  the 
beam  leaves  the  channel  the  tracking  force  drops  quickly  (see  fig.1).  Even  if  the  beam  is 
not  driven  out  of  the  channel  the  tracking  effect  can  be  masked  by  moderate  hose 
oscillations.  Ideally,  a  completely  hose  stable  beam  should  be  used,  but  for  a  demonstration 
of  the  tracking  effect  a  beam  which  is  stable  for  over  a  tracking  length  is  sufficient.  A  scale 
length  for  the  growth  of  the  hose  instability  is  the  betatron  wavelength  given  roughly  by: 


■J 


leff=  2J2 1 


Lb  =  2  n  lb 


(2) 


where  Ia  is  the  Alfv6n  current  (17f3y  kA),  |3  =  ‘Vc  is  the  normalized  electron  velocity  and 
Y=(1-{32)'^  is  the  relativistic  scaling  factor.  For  Pulserad  y=3.  Figure  3b  shows  the 
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calculated  betatron  wavelength  as  a  function  of  beam  radius.  For  the  Pulserad  beam 
Lb  »20-40cm. 

A  third  phenomena  which  can  negate  a  positive  tracking  force  is  electron  avalanche  in  the 
low  density  channel.  When  the  inductive  electric  field  produced  by  the  rising  current 
exceeds  the  avalanche  ionization  threshold  level  («  50kV/cm-atmosphere  for  air),  the 
collisionally  produced  plasma  electrons  gain  enough  energy  from  the  field  to  produce 
avalanche  ionization.  This  threshold  is  lowest  on  the  reduced  density  channel  axis.  If  the 
channel  avalanches,  the  additional  plasma  return  current  will  drive  the  beam  out  of  the 
channel.  A  simplified  expression  for  the  inductive  electric  field  near  the  head  of  the  beam 
is: 


Ez 


rb 

\  J 


kV/cm, 


(3) 


where  for  these  experiments  Tr(ns) » 10ns  is  the  risetime  of  leff(kA)  and  rw  is  the  wall  ra¬ 
dius,  (rw/rb  *25).  Figure  3c  shows  a  family  of  curves  of  the  axial  field  normalized  to  the 
avalanche  threshold.  At  high  currents  and  low  channel  density  the  threshold  is  ex¬ 
ceeded  which  should  result  in  repulsion  of  the  beam  from  the  channel. 


A  similar  effect  could  result  due  to  the  presence  of  free  electrons  in  the  channel  prior  to 
beam  injection.  The  LGED  method  used  to  produce  the  rarefied  channels  leaves  the 
channel  ionized  with  a  non-zero  conductivity  for  a  short  time  after  production.  If  the  DC 
conductivity  level  is  more  than  «  5x1 09  sec'1  at  beam  injection,  there  are  enough  free 
electrons  in  the  channel  for  the  beam  to  drive  sufficient  on-axis  return  current  to  expel  the 
beam  from  the  channel.[4]  The  LGED  channel  is  considerably  above  this  threshold  value 
when  first  created.  The  channel  cools  by  adiabatic  expansion  during  the  first  few  hundred 
microseconds  and  by  turbulent  convective  mixing  of  background  gas  into  the  hot  channel 
on  a  several  millisecond  timescale.  The  colder,  higher-density  gas  near  the  edges  moves 
toward  the  channel  axis  as  highly  asymmetric  three-dimensional  flutes  with  a  scale  length 
on  the  order  of  the  channel  radius.  Depending  on  the  total  energy  dumped  into  the  channel, 
2ms  or  more  must  elapse  before  beam  injection  to  prevent  the  beam  from  being  driven  out 
of  the  channel. 

Lastly,  higher-order  chemistry  processes  not  considered  in  the  theory  can  produce 
detracking  late  in  the  beam  pulse.  These  processes  include  electron-ion  recombination, 
electron  attachment,  and  Spitzer  collisions.  This  detracking  is  usually  offset,  however,  by 
longitudinal  magnetic  coupling.  Coupling  causes  the  beam  body  to  follow  the  head  and 
dominates  once  the  plasma  conductivity  is  large. [4]  Thus,  if  the  head  tracks  the  channel, 
the  body  and  tail  should  also.  The  higher-order  chemistry  processes  do,  however,  reduce 
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the  peak  tracking  force,  thereby  producing  tracking  distances  longer  than  predicted  by  the 
simple  analytic  theory.  Slinker[7]  and  Keeley[8]  have  conducted  detailed  numerical  studies 
of  the  effects  of  higher-order  chemistry,  including  water  vapor,  on  tracking. 

IV.  EXPERIMENTAL  TECHNIQUE 

Figure  4  shows  schematic  views  of  the  interaction  between  the  REB  and  a  displaced 
reduced-density  channel  in  the  experimental  setups  used  for  the  two  data  runs.  The 
chamber  was  «  90cm  inner  diameter  by  2.5m  long  and  made  of  lucite  lined  with  copper 
screen  to  facilitate  optical  photography  during  the  experiment.  The  electron  beam  from 
Pulserad  was  injected  into  the  experimental  chamber,  on  axis,  through  a  beam  current 
monitor.  Baseline  beam  propagation  data  was  obtained  by  injecting  the  beam  into  the 
chamber  with  no  channel  present.  In  the  first  experimental  run  a  9-1 2kA  beam  was  injected 
through  a  3.0mi I  titanium  (Ti)  anode  foil  into  ±1  cm  offset  channels.  The  net  current  centroid 
was  monitored  with  a  single  set  of  four  magnetic  probes  mounted  on  the  chamber  wall 
40cm  downstream  from  the  injection  point.  For  the  second  experimental  run  beam  currents 
were  decreased  to  5-1  OkA  at  injection  and  the  beam  radius  was  enlarged  to  about  2.1cm 
by  passing  it  through  an  additonal  1.5mil  Ti  foil.  The  reduced  current  and  larger  radius 
combined  to  increase  the  betatron  wavelength  (Eq.  2)  which  improved  the  beam  stability. 
The  net  current  centroid  position  was  measured  by  three  sets  of  magnetic  probes  20cm, 
60cm  and  110cm  downstream  of  the  injection  point.  The  rarefied  channel  was  offset  by 
±1 .4cm  and  by  ±2.6cm  from  the  chamber  axis,  or  placed  directly  along  the  chamber  axis. 
To  monitor  for  systematic  errors  in  the  beam  injection,  measurements  were  taken  in  groups 
of  four  successive  shots:  no  channel  present,  channel  displaced  up,  channel  displaced 
down,  channel  centered  on  chamber  axis.  Figure  5  shows  a  typical  beam  current  pulse 
measured  at  injection  and  the  net  current  20cm  downstream  for  a  no-channel  shot.  Figure 
6  shows  the  tracking  geometry  as  viewed  along  the  propagation  axis  for  an  offset  density 
channel  and  an  electron  beam  passing  along  the  edge  of  the  channel. 

V.  CHANNEL  PARAMETERS 

Extensive  measurements  were  made  on  LGED  channels.  Schlieren  photography  indicated 
that  the  channel  edge  radius  increased  from  below  3cm  to  about  4cm  in  the  time  between 
2ms  and  5ms  after  initiation.  Turbulent  mixing  with  the  surrounding  gas  appears  to  be  the 
primary  mechanism  for  decay  of  the  channels.  This  is  characteristic  of  using  an  aerosol 
and  a  laserto  designate  the  initial  breakdown  channel  in  the  LGED  process.  The  turbulence 
produces  a  region  on  the  sides  of  the  channel  with  a  mixture  of  high  and  low  density  gas. 
The  conductivity  generated  by  the  beam  passing  through  this  turbulent  region  should  be 
a  mixture  of  the  high  and  low  density  conductivities  in  series.  A  beam  produced  channel 
would  not  suffer  from  as  strong  turbulent  mixing.  The  channel  depth  infered  by  Stalder  et 
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al.[6]  from  early  time  conductivity  measurements  and  from  optical  interferometry  is  <,  Vio  at¬ 
mosphere.  The  depth  of  the  channel  and  the  diameter  are  dependent  on  the  bank  energy 
deposited  in  the  LGED  process. 

VI.  CHANNEL  TRACKING  EXPERIMENTAL  RESULTS 

Vla.--(First  run;  Summer  1987)The  9-12kA  Pulserad  beam  used  for  this  run  was  subject 
to  strong  hose  instability.  On  many  of  the  shots  the  beam  diverted  toward  the  wall  after 
propagating  to  about  Z=60cm.  The  magnetic  pickup  loops,  mounted  on  the  wall  at  Z=40cm, 
rw=45cm,  necessarily  averaged  over  the  beam  motion  for  some  distance  along  the  Z-axis. 
While  it  was  possible  to  determine  which  direction  the  beam  moved  if  it  hosed,  the 
magnitude  of  the  motion  was  an  average.  With  this  caveat  in  mind,  figure  7  presents  the 
net  current  centroid  measurements  for  this  run.  The  delay  time  between  channel  creation 
and  beam  injection  was  500-3500ps.  The  shot  data  has  been  normalized  to  the  direction 
of  the  channel  offset  such  that  motion  towards  the  channel  is  positive  and  motion  away 
from  the  channel  is  negative.  The  shots  with  no  channel  present  are  the  last  set  shown  in 
Fig.  7  and  are  designated  by  the  "NO-CH"  label.  At  each  delay  time  the  data  is  subdivided 
to  show  the  motion  in  the  beam  head  (1  Ons),  beam  body  (20ns)  and  the  beam  tail  (30ns). 
The  data  shows  the  beam  head  pulled  in  first,  the  the  body  and  tail  follow,  as  expected. 
The  error  bars  represent  one  standard  deviation  about  the  average  of  the  shots  taken  at 
that  particular  delay  time.  For  time  delays  <  2ms  the  beam  was  ejected  from  the  channel, 
probably  due  to  thermal  ionization  in  the  still  hot  channel.  After  a  delay  of  2.5ms  or  more, 
the  beam  moved  toward  the  channel  indicating  that  the  channel  had  cooled  enough  to 
permit  the  tracking  force  to  come  to  the  fore.  Measurements  of  the  bulk  channel  DC 
conductivity  indicated  that  it  dropped  below  the  theoretically  predicted  threshold  of 
5x1 09  sec*1  for  avalanche  ionization  at  that  time.  Since  the  no-channel  shots  were  not 
grossly  unstable  until  well  past  the  probe  position,  there  is  evidence  that  the  channel  was 
triggering  increased  hose  instability.  On  several  shots  the  beam  seemed  to  move  com¬ 
pletely  through  the  channel  on  its  way  to  the  wall,  consistent  with  numerical  simulations 
by  Taylor  et  al.[9].  Theoretical  predictions  for  the  tracking  length  of  a  >10kA  beam  with  a 
beam  radius  of  about  1.5cm  and  a  1cm  offset  channel  yield  lengths  less  than  100cm.  The 
data  from  this  high  curent  first  run  hint  at  a  short  tracking  length  but  the  hose  motion  of  the 
beam  made  it  impossible  to  extrapolate  a  valid  number  from  the  data.  Simulations  by  Taylor 
et  al.  agree  with  the  experimental  results.[9] 

Vlb.-(Second  run;  Summer  1988)The  5-1  OkA  Pulserad  beam  used  for  this  run  was  much 
more  stable  because  the  beam  current  was  reduced  and  the  beam  radius  increased  from 
the  first  run.  These  changes  also  increase  the  tracking  distance  Zt,  however,  as  was 
observed.  Three  sets  of  magnetic  probes  were  mounted  on  the  wall  (rw=45cm)  at  Z=20cm, 
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60cm  and  1 1 0cm  in  order  to  more  accurately  track  the  trajectory  of  the  beam.  It  was  found 
that  for  injected  currents  £  9kA  the  beam  was  hose  stable  for  the  entire  length  of  the 
propagation  range.  Figure  8  presents  for  the  no-channel  shots  the  average  beam  trajectory 
at  three  times  during  the  pulse.  The  error  bars  represent  one  standard  deviation  about  the 
average.  Since  the  channel  offsets  were  ±1.4cm  and  ±2.6cm  the  beam  motion  was 
acceptably  small.  As  with  the  first  run  the  beam  was  ejected  from  a  hot  channel  (delay  time 
£2ms).  Figure  9  illustrates  this  phenomenon  with  the  average  trajectory  for  beams  injected 
into  the  chamber  with  a  2.6cm  offset  channel  and  short  delay  times.  The  channel  size  and 
position  are  represented  by  the  shaded  region  in  the  figure. 

For  delay  times  £  3ms  (and  peak  beam  currents  still  <,  9kA)  the  trajectories  for  shots  with 
the  channel  along  the  chamber  axis  were  little  different  than  those  shots  without  a  channel 
present.  The  average  over  all  the  shots  in  the  set,  shown  in  figure  10,  indicates  that  the 
beam  trajectories  track  exactly  on  the  channel  axis  and  the  spread  of  the  trajectories  was 
not  much  larger  than  that  of  the  no  channel  shots.  The  small  difference  between  these  two 
data  sets  indicates  that  the  channel  did  not  induce  more  than  a  small  increase  in  the  hose 
instabliity  of  the  beam. 

The  next  two  figures  show  the  beam  motion  in  the  presence  of  oifset  channels,  again  for 
delay  times  £  3ms  and  injected  currents  <  9kA.  Figure  1 1  shows  the  average  trajectory  for 
the  ±1 .4cm  offset  channel  shots.  There  was  a  small  movement  of  the  beam  of  about  1cm 
toward  the  channel,  principally  at  the  peak  of  the  current  pulse  where  the  tracking  force  is 
strongest,  but  the  scatter  in  the  data  made  it  impossible  to  extrapolate  a  tracking  length 
from  this  data  set.  The  analytic  theory  predicts  the  tracking  force  should  increase  as  the 
offset  between  the  beam  and  channel  increases,  out  to  the  channel  edge,  but  the  tracking 
length  itself  increases  monotonically  as  the  offset  increases.  This  prediction  was  supported 
by  the  2.6cm  offset  channel  data  presented  in  figure  1 2.  The  average  beam  offset  at  peak 
current  was  more  than  2cm  towards  the  channel  center,  with  acceptably  small  scatter  in 
the  data,  at  the  Z=1 10cm  probe  position.  Extrapolating  the  average  trajectory  at  peak 
current  through  the  20cm,  60cm  and  1 1 0cm  data  points  to  the  channel  axis  yields  a  tracking 
length  of  Zt  - 1 60cm.  This  value  for  Z\  is  consistent  with  current  numerical  simulations.  The 
channels  are  about  Vio  th  density  at  center  shortly  after  being  created  but  they  fill  in 
asymmetrically  due  to  turbulence  as  the  channel  cools.  The  profiles  are  somewhat  ragged 
due  to  this  convective  mixing  of  cold  outside  gas  into  the  hot  channel  which  makes  for  a 
higher  average  density  and  an  increase  in  the  tracking  length  due  to  a  reduction  in  the 
attractive  force.  From  figure  1  a  reasonable  match  to  the  present  experiment  is  the 
calculation  of  the  interaction  between  a  5kA  beam  with  a  Bennett  radius  of  2.3cm  and  a 
2.6cm  offset  channel  with  a  half  density  radius  of  2.5cm  and  a  channel  density  about  0.5 
that  of  ambient  which  yields  a  calculated  tracking  length  of  150cm.  This  is  fairly  good 
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agreement  between  theory  and  experiment  considering  the  difficulty  in  determining  even 
an  approximate  profile  of  the  density  channel  for  use  in  the  calculations. 

VII.  DISCUSSION  AND  SUMMARY 

Both  experiments  demonstrated  clearly  the  existence  of  a  density  tracking  force,  as  well 
as  a  detracking  force  at  high  levels  of  cnannel  preionization.  Furthermore,  the  transition 
from  detracking  to  tracking  occurred  at  the  expected  channel  overheat  condition.[4] 
Although  the  tracking  length  and  overall  beam  motion  agreed  qualitatively  with  theory  and 
numerical  simulations, [9]  precise  quantitative  comparison  was  hampered  by  uncertainties 
in  the  experimental  parameters  and  by  complications  from  the  hose  instability. 

The  largest  uncertainties  are  the  channel  density  profile  and  the  beam  radius.  The  on-axis 
channel  density  was  *  0.1  atm  shortly  after  being  created,  but  then  gradually  rose  as  the 
channel  cooled  and  filled  in  asymmetrically  due  to  convective  mixing  with  surrounding  air. 
Although  the  complicated  spatial  structure  of  the  channel  makes  its  density  profile  difficult 
to  characterize  and  diagnose,  average  values  can  be  inferred  using  particle  conservation 
and  measured  values  of  the  (ragged)  channel  edge  radius.  The  beam  radius  is  similarly 
not  well  known  but  presumably  was  considerably  larger  in  the  second  set  of  experiments, 
as  indicated  by  the  longer  tracking  distances  and  more  stable  beam  behavior. 

Hose  instability  precluded  demonstration  of  improved  beam  propagation  in  the  reduced- 
density  channels.  We  hope  to  avoid  this  problem  in  future  experiments  using  the  Su- 
perlBEX  beam  by  conditioning  the  beam  prior  to  propagation.  Improved  diagnostics  will 
also  be  employed  to  provide  better  comparison  with  theory. 
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APPENDIX  I 


The  non-linear  form  of  the  tracking  force  equation[4]  can  be  evaluated  under  certain 
conditions.  Assume  the  total  electric  field  rapidly  approaches  the  axial  monopole  field  Ezo 
and  the  mobility  is  a  function  of  the  reduced  electric  field,  p  =  %©  =  .  The  angular 

part  of  the  force  equation  integral  evolves  to 


Here  the  form  of  the  mobility  is  assumed  to  be 

(E/n  V 

1-<Jln  —*  ( A2 ) 

lX°jl 

which  is  fit  to  published  values  of  the  mobility  of  air  using  q=0.2  and  xo=100  Townsends. 
The  calculation  is  closed  by  assuming  a  channel  density  profile  given  by 


ng(i')  =  floo 


V/acf  +  8  ip 
V/ac)2+  1 


(A3) 


h*  =  ambient  gas  density 

(O2  =  z2  +  y§  -  2  ryc  cos0 
Yc  =  channel  offset 
where  •  ac  =  channel  radius 
gP  =  channel  depth 

p  e  [0<p^  10] 

8  <  1 

Since  an  exact  channel  profile  could  not  be  determined  a  range  of  values  for  ac,  p  and  6 
were  chosen  for  a  parameter  study. 

The  partial  derivative  in  (A1 )  can  be  expanded  as  follows: 


_ 2 pyc  r  si nQ _  _ 2pyc  rsinQ _ 

/2  +  y|  +  8ac-2ycrcos0  z2  +  y§  +  ac-2ycr  cos© 


(A4) 


ll 


The  tracking  force  equation,  Ft  =  qE,  can  be  integrated  and  yields  an  asymptotic  tracking 
force  given  by: 


Ft  -  - 


2e  n 
ycc 


r  -*> 

'ne' 

r  r 

i 

,n9. 

4(1-6) 

Eza  +  [r4  +  2Z2  (6ai  -  y§)  +  (y§  +  ai)] 

. -[/  +  2^(4 -y§)  +  (y§+4)2]/2  J 


V& 


(A5) 


where  the  plasma  current  is  Ip  =  Jb  -  left  -  0.3  lb  at  the  peak  of  the  Pulserad  beam  pulse. 

ne  1 


Direct  production  causes 


to  be  self-similar  with  Jb.  Then  assuming  a  Bennett  beam 


with  radius  ab  one  can  integrate  (A5)  to  obtain: 
Ft  =  2kf^c][gi(1)"9l(5)]  Wher& 


51(5)  = 


(2ap  ac  yc  )* 


8ln(a)  + 


yc  -  (4  -  54  )2 


d  (8)  =  [  (y?+ 8a?  )2  +  24  (4  -  54 )  +  at] 

(d+84-4-4)4 


a  (8)  = 


d(4+64)  +  (yi+54)2  +  4(y§  -  §4 ) 


/c(8)  =  1  -  Q  In! 


*o 


M 


f  [In(8)  +  g2(1)-p2(5)] 


52(8)  = 


_  4+s4+4 


In  (a)  . 


(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 


The  tracking  distance,  Zt,  which  represents  the  propagation  length  necessary  for  a  beam 
to  be  pulled  to  the  channel  axis  is: 


Cl  _ 

ji 

Y  m(?  yc 

4  * 

2 

F< 

(A12) 


which  is  independent  of  yc  only  for  yc  ->  0. 
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Figure  1.  Pulserad  Simulation  of  the  tracking  force  and  the  tracking  length  for  Relative  Channel 
Density=0.5,  Channel  Radius=2.5cm,  Beam  Radius=2.3cm  andy=3. 
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Figure  2.  Approximate  tracking  length  using  nonlinear  theory  for  the  Pulserad  beam  with  n>=2.3cm  and 

rc=2.5cm 
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Figure  3.  Calculated  Pulserad  Beam  Propagation  and  Stability  Parameters. 
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Figure  4.  Schematic  diagram  of  the  two  propagation  chambers  used  in  the  two  experimental  runs. 
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Figure  6.  Pulserad  tracking  geometry  as  viewed  along  the  propagation  axis. 
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Figure  7.  '87  Channel  Tracking  Experiment  results  for  net  current  centroid  motion  versus  delay  time  after 

channel  creation. 
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Figure  8.  Beam  Trajectory  for  the  Pulserad  Only,  No-Channel  Shots  with  Beam  Current  <,9kA. 
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Figure  9.  Beam  Trajectory  for  the  2.6cm  offset  Channel  Shots  where  the  delay  time  between  channel 

creation  and  beam  injection  was  <,2ms. 
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Figure  1 1.  Beam  Trajectory  for  the  1.4cm  offset  Channel  Shots  with  the  beam  current  <,9kA  and  the  delay 

time  z3ms. 
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10529B  Braddock  Road 
Fairfax,  VA  22032 
Attn:  Dr.  F.  M.  Mako 


GA  Technologies,  Inc. 

P.  O.  Box  85608 
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San  Diego,  CA  93138 
Attn:  Dr.  Vincent  Chen 

Dr.  Hiroyuki  Ikez 

General  Dynamics  Corporation 
1745  Jefferson  Davis  Highway 
Suite  1000 

Arlington,  VA  22202 
Attn:  Dr.  Daniel  W.  Roth 

General  Dynamics  Corporation 

Pomona  Division 

1675  W.  Mission  Blvd. 

P.  O.  Box  2507 
Pomona,  CA  92769-2507 
Attn:  Dr.  Ken  W.  Hawko 

Mr.  C.  L.  Featherstone 

Grumman  Corporation 
Grumman  Aerospace  Research  Ctr . 
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Dr.  Francis  Tsang 

27 


Institute  for  Defense  Analyses 
1801  N.  Beauregard  Street 
Alexandria,  VA  22311 
Attn:  Dr.  Deborah  Levin 
Ms.  M.  Smith 

IRT  Corporation 
3030  Callan  Road 
San  Diego,  CA  9 2121 
Attn:  Dr.  David  Phelps 

JAYCOR 

39650  Libery  Street,  Suite  320 
Freemont,  CA  94538 
Attn:  Dr.  Kendal  Casey 

Joint  Institute  for  Laboratory 
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National  Bureau  of  Standards  and 
University  of  Colorado 
Boulder,  CO  80309 
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Dr.  George  J.  Caporaso 

Dr.  Donald  Prosnitz 

Dr.  Y.  P.  Chong 
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3251  Hanover  St. 

Bldg.  205,  Dept  92-20 
Palo  Alto,  CA  94304 
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Attn :  Dr .  L .  Thode 
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William  Barletta 

Dr. 
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Dr.  Michael  Mostrom 
Dr.  Dale  Welch 

Mission  Research  Corporation 
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Naval  Surface  Warfare  Center 
White  Oak  Laboratory 
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Dr.  H.  C.  Chen 
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Office  of  Naval  Research 
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Attn:  Dr.  C.  W.  Roberson 
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Department  of  the  Navy 
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Attn:  Dr.  Jose  F.  Pina 
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Woburn,  MA  01801 
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2700  Merced  Street 
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Dr.  James  Benford 
Dr.  George  B.  Frazier 
Mr.  Ralph  Genuario 

Princeton  University 
Plasma  Physics  Laboratory 
Princeton,  NJ  08540 
Attn:  Dr.  Francis  Perkins,  Jr. 

Pulse  Sciences,  Inc. 
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Attn:  Dr.  Sidney  Putnam 
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Alexandria,  VA  22314 
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Dr.  D.  Chernin 
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2755  Northup  Way 
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Attn:  COL  Thomas  Meyer  (DEWO) 
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University  of  California 
Physics  Department 
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Los  Angeles,  CA  90024 
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University  of  Colorado 
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Dr.  W.  Destler 
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Yale  University 
Mason  Laboratory 
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Attn:  Dr.  Ira  Bernstein 
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University  of  Illinois  at  Chicago 

Dept,  of  Physics 
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Attn:  Dr.  Charles  K.  Rhodes 
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